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Non-linear saturation of g-modes in proto-neutron stars: 
quieting the acoustic engine 

Nevin N. Weinberg, 1 * Eliot Quataert 1 

1 Astronomy Department and Theoretical Astrophysics Center, University of California, Berkeley, 601 Campbell Hall, Berkeley CA, 94720 



00 

o 

O 

(N 

C 



43 

Oh 
6 



> 

(N 

in 

O 

oo 
o 



X 



Accepted. Received; in original form 



ABSTRACT 

According to Burrows et al.'s acoustic mechanism for core-collapse supernova 
explosions, the primary, 1 = 1, g-mode in the core of the proto-neutron star is excited to 
an energy of ~ 10 50 ergs and damps by the emission of sound waves. Here we calculate 
the damping of the primary mode by the parametric instability, i.e., by nonlinear, 3- 
modc coupling between the low-order primary mode and pairs of high-order g-modes. 
We show that the primary mode is strongly coupled to highly resonant, neutrino 
damped pairs with n > 10; such short wavelength interactions cannot be resolved in 
the simulations. We find that the parametric instability saturates the primary mode 
energy at ~ 10 48 ergs, well below the energy needed to drive an explosion. We therefore 
conclude that acoustic power is unlikely to be energetically significant in core-collapse 
supernova explosions. 

Key words: hydrodynamics — instabilities; stars: neutron — supernovae: general 



1 INTRODUCTION 

Of the ~ 10 53 ergs of gravitational binding energy released 
in a stellar core-collapse, ~ 1% is used to power the su- 
pernova explosion and the remainder is nearly all lost to 
neutrinos that escape to infinity. According to the neutrino 
mechanism, the supernova is driven by the small fraction 
of neutrinos that deposit energy in the gain region just be- 
low the stalled accretion shock. Over 20 years of study has 
demonstrated that the viability of the neutrino mechanism 
hing es on the multidim ensional nature of core-collapse (see 
e.g., I Janka et al.l 1200*71 for a review), wit h some recent 2D 
simulations on the borderline of success l|Buras et alj|2006l ; 
iMarek fc Jankall2007l ). 

Recently, Burrows et al. (2006, 2007) proposed a new 
explosion mechanism that instead relies on acoustic power 
generated in the core of the proto-neutron star (PNS). Us- 
ing 2D, axisymmetric, multigroup, radiation-hydrodynamic 
simulations, Burrows et al. find that after the stalled shock 
becomes unsta ble to the standing accretion shock instabil- 
ity (S ASI; see iBlondin fc Mezzacappal 120061 ; iFoglizzo et~ai1 
120071 ). anisotropic accretion streams excite and maintain vig- 
orous g-mode oscillations in the PNS core. Sound pulses ra- 
diated from the core steepen into shock waves and deposit 
sufficient energy and momentum in the outer mantle to drive 
an asymmetric explosion ~ 1 s after bounce. 

The primary mode of oscillation seen by Burrows et 
al. is a low-order core g-mode with spherical degree / = 1 
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(frequency « 300 Hz). For an IIM0 progenitor, they find 
that the pulsational energy in the core is maintained at 
~ 10 50 ergs for several hundred milliseconds before the ex- 
plosion and that the dominant source of damping is the emis- 
sion of sound waves near the PNS surface. For more massive 
progenitor models, the core pulsations achieve energies near 
~ 10 51 ergs. 

As the primary mode is driven to ever larger amplitude, 
nonlinear coupling to other PNS modes can become an im- 
portant source of damping. Here we calculate the leading- 
order nonlinear coupling between the primary mode and 
high-order core g-modes (the parametric instability). The 
latter rapidly damp due to neutrino diffusion. Burrows et 
al.'s resolution in the PNS core is w 0.5 km and they can- 
not adequately resolve the coupling of the primary to higher 
order modes; an analytic calculation is thus required. 

We begin in § [2] with a brief review of the parametric 
instability and outline the calculation of the saturation en- 
ergy. We then describe our PNS model and its eigenmodes 
(§ I2.1[1 . their coupling strength (§ I2.2[1 . and their driving and 
damping rates (§§ 12.31 T2.4p . In §[3] we show results and dis- 
cuss their implications. A more detailed presentation of our 
methods and results will be presented in a separate paper. 



2 SATURATION BY NONLINEAR MODE 
COUPLING 

In the particular context of three-mode coupling con- 
sidered here, if the amplitude of the I = 1 "parent" 
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mode exceeds the parametric instability threshold, it be- 
gins to tra nsfer energy to tw o lower-frequency "daugh- 



ter" modes llDziembowskil 1 19821 ; iKumar fc Goodman] 19961 : 



IWu fc Goldreichll200ll . hereafter WG01: I Arras et alj(200 j ) 
The daughter mode amplitudes then rise exponentially, and 
the system has an equilibrium in which the parent mode 
energy is 



Ei „q — 



7273 



18k 2 LU2U>3 



1 + 



25uj 



72 + 73 — 71 



(1) 



Here Wj is the eigenfrequency of mode i (an index of 1 
refers to the parent mode and 2 and 3 to the two daugh- 
ter modes) , Slo = uji — L02 — u)z is the detuning frequency, 
72,3 are the daughter damping rates, 71 is the parent driving 
rate, and k is the nonlinear coupling coefficient The daugh- 
ter modes' equilibrium energy is i?2,3,eq/-Ea,eq = Q2, 371/wi, 
where Q2,3 = 1^2,3/72,3 is the daughter mode quality factor. 
As we describe in §[3] the parent's saturation energy _E sa t is 
set by the daughter pairs that minimize -Ei, e q- 

Accretion onto the PNS pumps energy into the primary 
mode at a rate 71 = E\/E\, where according to Burrows et 
al. Ei ~ 10 50 - 10 51 ergs s . Unlike in many other applica - 
tions of the parametric instability (e.g., lArras et al.ll2003h . 
the driving here is nonlinear because Ei, rather than 71, 
is taken to be constant. Equation (jXJ is thus a cubic equa- 
tion in Ei. For \5uj\/lui < 0.01, it has three real roots if 
£1/(72+73) > [lS^QaQg]" 1 and one real root otherwise. 
In the three root case, the largest root -Bi, eq ~ £'1/(72 + 73) 
is the stable solution (as revealed by both an analytic sta- 
bility calculation and direct numerical solution of the time- 
dependent amplitude equations for the 3- mode system). 
In the one root case, the solution for small detuning is 

£i,eq « [lS/^QaQsr 1 . 

Before describing the calculation of k and 7, we first 
give an estimate of E aa ,t- In our units, the total energy of a 
mode is normalized to unity and the coupling coefficient for 
core g-modes is of order k ~ l/V-Ebind = (GM 2 /R)~ 1/2 « 
(10 53 ergs) _1/2 , the inverse square-root of the binding en- 
ergy of a PNS of mass M and radius R. The strength of 
the 3-mode coupling is sensitive to the structure of the 
modes, particularly the turning points and the location of 
nodes. While for most daughter pairs the coupling is weak 
(k <C l/V-Ebind), we show that some pairs couple strongly to 
the primary, with k » lO/V-Ebmd- Of the strongly coupled 
pairs a few with n/l ~ 5 have damping rates 72,3 w 100 Hz 
and |5w|/wi ~ 0.01. For these fiducial daughter modes, if 
Ei < 10 50 ergs s _1 , the primary saturates at an energy of 



10 4 



W^/Rao/Ml 



(9m. 
\ 10 



ergs, 



(2) 



(.R30 = R/3Q km and Mia = M/1.4M©) while for larger Ei 
the three root solution to equation JTJ applies and E &a x ~ 
5 x 10 48 (i;i/10 51 ergs s" 1 ) ergs. 

After the parent crosses the parametric threshold, 
the energy growth rate of the daughter modes is (72,3 w 



[I8E1K LU2UJZ — {810) 



211/2 



72,3. Thus, if the parent over- 



shoots its equilibrium by just a factor of two (Ei ~ 2_Ei, c 



We follow the normalization conventions of WG01. 




5 10 
radius (km) 

Figure 1. Profiles of the radial displacement g r {top panel; nor- 
malized such that the total pulsational energy is one), daughter 
damping rates (middle panel; 7(< r) = / (d-yfdr)dr), and cou- 

pling coefficient (bottom panel; k(< r) = J (dn/dr)dr, in units 

of 10~ 26 ergs -1 / 2 ) for 3-mode coupling of the primary mode to 
a pair of modes with (l,n) = (2,14) and (3,9). Table [T] lists ad- 
ditional properties of this triplet. 



a ~ 100 Hz and the daughters undergo many tens of e- 
foldings of growth in just a small fraction of a second. Al- 
though the initial daughter mode amplitudes are uncertain, 
we expect PNS convection (and perhaps also the asymmetric 
accretion streams) to excite them to energies S> 10 40 ergs. 
Equilibrium [equations (JTJ & (|2J] is therefore established in 
< 0.1 sec. At the energies found in Burrows et al.'s simula- 
tions (Ei ~ 10 50 ergs), the daughters would grow in ~ 1 ms 
and very quickly drain energy out of the primary. 



2.1 Proto-Neutron Star Model and Eigenmodes 

The PNS models used in our study were provided by A. 
Burrows. They are angle- averaged radial slices of density 
(p), temperature (T), electron fraction (Y e ), and pressure 
(p), at variou s times during their core-collapse simulation 



of the 15M 



o 



iBurrows et al 



Wooslev et al. (|2002T ) progenitor model (see 



20071 ). We show results for the PNS at t 
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Table 1. Daughter pairs with low parametric thresholds 



(h,n 2 ) ■ 


(l:i,n 3 ) 


\Su\/u>i 


72,73 




-Ei, cq (ergs) a 


(8,39) : 


(9, 33) 


0.5 


132,69 


1.2 


5.1 X 10 4S 


(3,21) : 


(4, 13) 


2.5 


740, 25 


1.3 


7.2 x 10 48 


(1,16) 


: (2,4) 


16.2 


548, 12 


0.9 


8.0 x 10 48 


(2, 14) : 


(3,9) b 


5.7 


241, 250 


1.8 


1.3 x 10 49 


Note. — |<5u|/u;i in 


10~ 3 ; 72,73 


in Hz; |re 


in 10" 


-26 ergs -1 / 2 . 



a Values of -Ei, e q are for E\ = 10 51 ergs s 1 . 
b See Fig. [T]for £, r (r), y(r), and n(r) for this triplet. 



500 ms after core bounce, although the value of E sat is nearly 
constant in timeQ 

We interpolate over the model and compute the 
adiabatic oscillation frequencies and eigenfunctions us- 
ing the non-relativistic Aarhu s adi abatic oscillation 
package (IChristensen-Dalsgaardl 120071 ). We use the 
lLattimer fc Swestv 1 199lf) equation of state (EOS) to cal- 
culate the Brunt- Vaisala (buoyancy) frequency N, chemical 
potentials p, etc., needed in the eigenfunction and damping 
calculations. Since the modes we consider all have damping 
rates 7 « w, the adiabatic approximation is appropriate. 
The eigenfunctions and coupling coefficient are computed 
in the Cowling approximation (i.e., we ignore perturbations 
to the gravitational potential), which is appropriate for 
the high-order daughter modes under consideration and 
should also be reasonable for the parent mode. We assume 
that the Lagrangian pressure perturbation vanishes at 
the outer boundary, which we place at p — 10 9 g cm~ 3 
(corresponding to a radius r- ou tcr — 80 km). Our results are 
not sensitive to the location of the outer boundary since the 
nonlinear coupling is dominated by the core (r < 15 km) 
and the mode damping is dominated by the core and the 
region near the neutrinosphere (r w 40 km). 

Motivated by Burrows et al.'s results, we assume that 
the primary (i.e., parent) mode has a spherical degree h — 1 
and radial order ni = 2. We find that it is a core mode 
with frequency Wi = 2-7T x 297 Hz, in good agreement with 
Burrows et al. . Although we use the Lattimer-Swesty EOS 
instead of the IShen et al.l (|l998l ) EOS used by Burrows et 
al., this agreement is reasonable given that we use Bur- 
rows et al.'s {p, T, Y e ,p} profiles; i.e., the Brunt- Vaisala fre- 
quency of the two EOS's — which largely determines u> — 
differ by ST 1 /2F 1 < 10%, half t he fractional difference o f 
the EOSs' adiabatic index Ti (see lRosswog fc Da vics 2002). 
Although the I = 1, n = 1 mode has a similar frequency 
(324 Hz), it is a surface mode that lies almost entirely above 
the PNS convection zone (r > 20 km; similar to Fig. 4 of 
lYoshida et all 120071 ). We find that it has a neutrino damp- 
ing rate of ~ 10 3 Hz. It therefore seems inconsistent with 
Burrows et al.'s primary mode and we instead focus on the 
n = 2 mode0 



We consider the 3-mode coupling of the primary to all 
daughter modes with 1 ^ I ^ 10 and frequencies in the 
range u)/uj\ ~ [1/4,3/4] (corresponding to 3 ^ n ^ 18 for 



I = 1 and 30 ^ n ^ 64 for I = 10). Selection rules enforce 
conservation of angular momentum ( | Za — Z3 1 ^ li ^ I2 + li 
with h + I2 + h even; see WG01), leaving ~ 6000 triplets 
with nonvanishing k. The top panel of Figure [T] shows the ra- 
dial displacement £r(r) of the primary mode and a strongly- 
coupled, resonant daughter pair. 



2.2 Coupling Coefficients 



The three-mode 
modes under the 



( Kumar fc Goldrcic 



WG01; Schenket 



ich 
al. 



coupling coefficient for adiabatic 
Cowling approximation is given by 



.198' 
120021 ) 



lg approximati on is gi v en by 
1 iKumar fc Goodman! Il996l ; 



d 3 xp{ 



(Ti - l) 2 + 



9r a 



dlnp 



(v • o 3 



+3 (r 1 -i)(v.0£4 + 2 £4^ 



+g / ^xpeee^jk, 



(3) 



where p, p, and <f) are the unperturbed pressure, density, and 
gravitational potential, Ti is the adiabatic index, and £ is 
the Lagrangian displacement vector. The semicolon denotes 
covariant derivative and the integration is over the volume of 
the star. As in appendix Al of WG01, we carry out the an- 
gular integrations analytically and use integration by parts 
to put this expression in a form more suitable for numerical 
computation|j The bottom panel of Figure [1] shows da/dlnr 
and k(< r) — f^(dn/dr) dr for the coupling of the primary 
mode to a strongly coupled, resonant daughter pair. 

The strength of the coupling depends on the triple prod- 
uct of the covariant derivative of each mode's eigenfunction 
and is very sensitive to the location of nodes and turning 
points (i.e., where N = ui). We find that even after we carry 
out the angular integrations analytically and use integration 
by parts, the modes need to be resolved to < 0.1 km in or- 
der to accurately calculate k. If we instead calculate ft by 
directly integrating equation @ over a cartesian grid, a res- 
olution of < 0.01 km is needed. It is difficult to infer from 
these analytic techniques what resolution a hydrodynamic 
simulation must have in order to accurately capture the non- 
linear interaction; we nonetheless believe that the resolution 
of current simulations (> 0.5 km) is not sufficient. 

If daughter modes parametrically excite granddaugh- 
ter modes before reaching equilibrium with the primary, 
they may saturate at e nergies below _E2,3,oq (see e.g., 
IKumar fc Goodman! [i"99rj| ; WG01). By searching the (l,n) 
parameter space for granddaughter modes that can couple 
to daughter modes, we find that daughters with I < 4 do 
not excite granddaughter modes at their equilibrium energy 
^2,3, cq ; for higher I daughter modes there generally are pairs 
of granddaughters that can be parametrically excited. Ac- 
counting for such coupling requires a full mode network cal- 
culation which we leave to a future paper. As can be inferred 



2 We also obtain E sa t ~ 10 48 ergs for a PNS model taken from a 
ID core collapse s imulation of an HM^ progenitor (provided by 
T. Thompson; see iThompson et al.ll2005h . 

3 We find that the n = 1 mode also has -E S at ~ 10 48 ergs. 



4 There are errors in terms 5, 7, and 9 of WGOl's final expression 
(their eq. [A15]; the errors are also present in eq. [A14]). We 
correct these and also include the unperturbed gravity term (last 
term in our eq. [3]), which they neglect. 
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from Table [T] however, when we limit our calculation to 
daughters with I < 4, we still find that E aa .t ~ 5x 10 48 ergs. 



2.3 Primary Driving 



Burrows et al. find that accretion pumps energy into the 
primary mode at a rate E\ ~ 10 50 — 10 51 ergs s _1 . In solv- 
ing the time-dependent coupling between the parent and 
daughters, we thus assume that the parent is driven with 
a constant energy injection rate of E\ = 10 50 ergs s" 1 or 
Ei — 10 51 ergs s _1 . A constant Ei corresponds to a nonlin- 
ear driving rate 71 = E1/E1 in equation ((TJ. 

Since the accretion streams that excite the primary are 
generated by the nonlinear SASI oscillation, an alternative 
modeling approach is to assume a constant linear driving 
rate 71 equal to the S ASI oscillation frequency. According to 
iBurrows et al.l l|2007l ). the frequency of the dominant, time- 
averaged, nonlinear SASI oscillation ranges from 30 Hz (for 
the 11.2M progenitor model) to 80 Hz (the 25M Q model). 
As can be inferred from equation (JTJ, such a linear driving 
model (e.g., with 71 = 50 Hz) yields very similar results to 
the nonlinear driving model. 



2.4 Damping Rates 

In the neutrino-degenerate PNS core, the energy damping 
rate of the daughter modes, 72,3, is dominated by neutrino 
damping due to conduction and particle diffusion (see e.g., 
Ivan den Horn fc van Weert| [l984)Fl The modes are evanes- 
cent in the PNS convection zone (13 km < r < 20 km) and 
daughter modes with wavelengths A much shorter than the 
size of the convection zone have negligible amplitudes out- 
side of the core. By contrast, longer wavelength daughter 
modes tunnel through the convection zone and have a non- 
negligible amplitude near the neutrinosphere at r v ~ 40 km. 
Such modes experience significant additional damping due 
to optically thin emission in the region rthin <r <r„ where 
the neutrino mean free path {d v ) > A (rthm is the radius 
where (d v ) = A). Modes with A ~ r v — r t hi n ~ few km 
experience significant optically thin damping at these radii, 
while longer wavelength modes do not. 

We account for damping in the core and near the neu- 
trinosphere as follows. We assume that the damping occurs 
in the quasi-adiabatic limit such that 7 = ~{E)/E to t where 
(E) is the time-averaged work integral expression for the 
energy lost and Etot is the total pulsational energy. We cal- 
culate the Rosseland mean free path of the ne utrinos (d„) 
using cross-sections given in lReddv et al.l (|l998l ). We assume 
7 = 7thick + Tthin, where 



7thick 



2-Eto 



-V • N v ) (4) 



accounts for optically thick damping in the core and 



7thin 



2E 



1 I™ 3 8T 

d xp—Se 



(5) 



r thi, 



5 Turbulent viscosity in the convection zone increases the damp- 
ing rates by at most as 30% and we therefore ignore it here. 




0.2 0.3 
time (s) 
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Figure 2. Energy of the driven primary (I = l,n = 2) g-mode in 
a PNS as a function of time for coupling between the primary and 
all g-modes with 1 ?C I ^ 10 and 1/4 < lj/oji < 3/4 (not including 
daughter only or granddaughter couplings). The primary mode is 
assumed to be driven with a constant energy injection rate of 



Ei = 10 50 ergs s 1 (bottom curve) and Ei 



10 s 



(top 



curve), corresponding to the power supplied by accretion seen 
in the simulations of Burrows et al. The initial energy of the 
daughters is 10 43 ergs. 



accounts for optically thin damping near the neutrinosphere. 
In equations Q & |[5j, 8 denotes a Lagrangian perturba- 
tion, N„ = -a T VT - a v Vr] and F v = -fc T VT- b v V77 are 
the neutrino number flux and energ y flux with transpor t 
coefficients {ax, a v , 6t, b v } (see e.g.. iMiralles et all feoOO), 
77 — p v jkT is the neutrino degeneracy parameter, and pt v 
is the neutrino chemical potential. The neutrino emissivity 
e is due to charge-current processes and the production of 
electron, muon, and tau neutrino pairs, which we calculate 
using the formulae in Thompson et al. [2001; equations (20) 
& (29)] and Thompson et al. [2000; equation (53)]. 

The middle panel of Figure [1] shows d r y/dlnr and 
7(< r) = f (cfaf/dr)dr in the core for a strongly-coupled, 
intermediate-order, resonant daughter pair. The majority 
(w 95%) of the damping for this pair occurs outside the 
core, in the optically thin region (see Table [TJ. By contrast, 
we find that the damping rate of the highest-order modes 
(n > 25) is dominated by the core and is to a good approx- 
imation, 7 w 7thick ~ 0.1 (n 2 + i 2 ) Hz; this is consistent 
with a PNS with a w 10 s Kelvin-Helmholtz cooling time. 



3 RESULTS AND DISCUSSION 

To accurately determine the saturation energy of the (I = 
1,71 = 2) parent mode, we solve the time-dependent am- 
plitude equations that describe the coupling of the parent 
mode to the set of high-order g-mode pairs described in ij l2.ll 
Figure [5] shows the resulting evolution of the parent mode 
energy for E\ = 10 50 ergs s _1 and Ei = 10 51 ergs s _1 . 

The parent mode energy saturates at E sa x ~ 1.5 x 



© 2008 RAS, MNRAS 000, [Tj{6] 



Saturation of g-modes in proto-neutron stars 5 




Figure 3. Equilibrium energy of the primary g-mode as a func- 
tion of detuning frequency \8u}\/u>i for individual 3- mode coupling 
between the primary and each g-mode pair with 1 ^ I ^ 10 and 
1/4 < u>/u>i < 3/4. We show results for E\ = 10 50 ergs s _1 (as- 
terisks) and Ei = 10 51 ergs s — 1 (circles). 



10 47 ergs for Ei = 10 50 ergs s" 1 and E S3 . t es 2 x 10 48 ergs 
for Ei — 10 51 ergs s . At these energies, the parent and 
daughter modes are only mildly nonlinear, with Sp/p < 0.1; 
the nonlinear interactions are therefore well-described by 
the leading-order nonlinear terms accounted for here. Al- 
though the simultaneous interaction between the parent and 
all daughter pairs is clearly more intricate than individual 
3- mode coupling, Figure [3] shows that the latter provides a 
good estimate of the parent's saturation energy. Note that 
there are dozens of daughter pairs with -Ei, cq ~ 10 48 ergs 
(5 x 10 48 ergs) for Ei = 10 50 ergs s" 1 (10 51 ergs s" 1 ). Pa- 
rameter values for a few of the low-_Ei je q daughter pairs are 
given in Table [T] 

Our calculations imply that most of the energy supplied 
to the parent core g-mode by accretion onto the PNS will 
be thermalized as heat via neutrino damping of higher order 
(n > 10), parametrically resonant daughter pairs; such short 
wavelength interactions cannot be resolved in current core 
collapse simulations. For a given saturation energy E sa t, the 
acoustic power radiated by the PNS is ~ 7 ac -E sa t, where 7 ac 
is the acoustic damping rate of the parent, which Burrows et 
al. find is 7 ac ~ 10 Hz. For simplicity we assume that 7 ac is 
the same at the low saturation energies found here[f| in which 
case the acoustic power radiated is ~ 10 48 — 2 x 10 49 ergs s _1 
for the values of _E sat in Figure [2] This corresponds to a few 
percent of the power supplied to the primary core g-mode. 
This level of energy injection into the surrounding star is 
insufficient to drive a core-collapse explosion. In order for 
acoustic power to generate a successful ~ 10 51 erg explosion, 
we find that > 1 - 3 x 10 52 ergs of accretion energy must 



6 This need not be the case if the coupling between the parent 
and outgoing sound waves is nonlinear. 



be supplied to the primary mode. This is ~ 10 — 100 times 
more energy than is seen in Burrows et al.'s simulations. 

The damping rates 72,3 and coupling coefficients k are 
subject to uncertainties that could influence the results pre- 
sented here. In particular, determining the exact value of 
7thin [equation ([5])] will require a more careful treatment of 
neutrino transport near the neutrinosphere. Uncertainties 
in the equation of state and the neglect of general relativ- 
ity and rotation introduce additional uncertainties into 72,3 
and k. To examine the effect of these uncertainties on the 
primary's saturation energy -E sat , we have recalculated the 
time-dependent amplitude equations with all of our fiducial 
k's and 7thin's multiplied by constant artificial factors. For 
Ei = 10 51 ergs s , the saturation energy is determined by 
daughters confined to the core and E sa _ t ~ £'i/7thick (the 
three root solution in § |2J. We thus find that even if we 
multiply/divide our k's and 7thin's by factors of ~ 3, i5 sa t 
only changes by 20 — 30%. Since the diffusive approximation 
is well-satisfied in the inner ~ 10 km, and since our fidu- 
cial 7thick's [equation Q] are consistent with a PNS with a 
Kelvin-Helmholtz cooling time of tun ~ 10 s, we do not be- 
lieve that 7thick is uncertain at the factor of ~ 3 level. In par- 
ticular, decreasing 7thick (and thus increasing -E sa t) by more 
than a factor of 3 would imply £kh > 30 s , inconsistent wit h 
detailed calculations of PNS cooling (e.g.. |Pons et al.lll999l ). 
For Ei ~ 10 50 ergs s _1 , daughters confined to the core again 
determine the saturation energy, but £ sat ~ [18k 2 Q 2 Q3}~ 1 
(the one root solution in § [2J). As before, because the core 
damping is in the diffusive limit, uncertainties in the rel- 
evant daughter quality factors are unlikely to be large. In 
addition, even a factor of 3 decrease in k only increases -E sa t 
to ~ 10 48 ergs. We thus believe that the general conclusion 
of this Letter is reasonably robust to uncertainties in PNS 
microphysics. In the future, a more detailed mode network 
will be utilized to further assess this conclusion. 

Several additional consequences of our results should be 
noted. The increase in neutrino flux due to the thermalized 
pulsational energy is negligible compared to the total flux 
diffusing out of the PNS. It is therefore unlikely to assist the 
neutrino mechanism. Our results also imply that g-modes in 
PNSs cannot be strong sources of gravitational wave radi- 
ation. Lastly, future core-collapse simulations must account 
for sub-grid mode coupling in order to suppress the artifi- 
cial growth of low-order core oscillations in cases where such 
growth influences the simulation results. 
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